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Motivated by recent nuclear magnetic resonance experiments on ZnCu3(OH)6Cl2, we present an exact- 
diagonalization study of the combined effects of non-magnetic impurities and Dzyaloshinskii-Moriya (DM) 
interactions in the s — 1/2 kagome antiferromagnet. The local response to an applied field and correlation- 
matrix data reveal that the dimer freezing which occurs around each impurity for D — persists at least up 
to D/J ~ 0.06, where J and D denote respectively the exchange and DM interaction energies. The phase 
transition to the (Q — 0) semiclassical, 120° state favored at large D takes place at D/J ~ 0.1. However, 
the dimers next to the impurity sites remain strong up to values D J, far above this critical point, and thus 
do not participate fully in the ordered state. We discuss the implications of our results for experiments on 
ZnCu3(OH)6Cl2. 

PACS numbers: 75.10.Jm, 75.3().Gw, 75.30.Kz, 75.40.Cx 



I. INTRODUCTION 

The s = 1/2 antiferromagnetic (AFM) Heisenberg model 
on the two-dimensional kagome lattice is one of the sim- 
plest models in frustrated quantum magnetism, but displays 
some of the most complex behavior known. ' As such, it has 
for decades maintained its position at the forefront in the 
search for novel quantum-mechanical phases of matter, such 
as the resonating-valence-bond (RVB) spin-liquid state pro- 
posed initially by Anderson.^ The nature of the ground state 
of the s = 1/2 kagome AFM has still not been fully estab- 
lished: the variety of competing phases proposed in the lit- 
erature includes valence-bond crystals (VBCs),^ gapped spin 
liquidsr^ and gapless critical phases. Under these circum- 
stances, the discovery that the kagome AFM is extremely sen- 
sitive even to the smallest perturbations, such as the presence 
of anisotropies'''"' or of non-magnetic impuritiesr^ii^ is com- 
pletely consistent. 

In this context, the recent discovery'^ of the mineral her- 
bertsmithite, ZnCu3(OH)6Cl2, has attracted very considerable 
attention, because it represents a structurally perfect realiza- 
tion of the s = 1/2 kagome AFM. This material is composed 
of Cu^+ ions (s — 1/2) arranged on kagome planes separated 
by triangular layers of non-magnetic Zn^+ ions. Despite the 
large AFM exchange, J ~ 170-190 K in this materialr^ 
in experiment there is no evidence of long-ranged magnetic 
order or even of spin freezing at any temperatures down to 
50 mKii^iii There is in addition no sign of a spin gap in the 
excitation spectrum. '^''^ 

However, it has also been established that in 
ZnCu3(OH)6Cl2 there is a significant (5-10%) intersite 
exchange of Cu^+ ions with the Zn^+ ions intended to 
separate the kagome planes. The displaced Cu^+ ions are 
thought to account for the large Curie tails observed in powder 
susceptibility measurements'^'^' and for the field-dependent, 
Schottky-like anomaly in the specific heat.^" At the same 
time, the Zn^+ ions displaced into the kagome planes play 
the role of non-magnetic vacancies, which are known" '^ 



to modify the ground state properties in a nontrivial way. 
However, these impurities can also be considered as a probe 
of kagome physics, and in this respect the recent '^O nuclear 
magnetic resonance (NMR) experiments reported by Olariu 
et al}^ offer extensive insight. While each O^^ ion is coupled 
predominantly to two neighboring Cu^+ ions in the kagome 
planes, the NMR spectra revealed a broad distribution of local 
susceptibilities. The relative 1:2 ratio between the line shifts 
of the two leading features in the '^O spectrum was explained 
on the basis of distinguishing between two groups of O sites 
in the doped kagome planes, those ["Defect" (D)] probing 
the magnetic polarization directly next to an impurity site, 
and hence coupled only to one Cu spin, and those ["Main" 
(M)] reflecting the polarizations of all other sites. This 
interpretation was also consistent with the relative intensities 
of the two lines based on the expected in-plane dopant 
concentration of approximately 5%. Another major finding 
of these experiments was the observation that the (M) and 
(D) line shifts approach a finite value as T ^ 0, suggesting a 
non-singlet ground state without a gap. This result can also 
be inferred from the non-activated behavior of the nuclear 
spin-lattice relaxation rate, 1 /Ti, at low temperatures. ' - 

While the true nature of the non-magnetic ground state 
of the pure, s — 1/2 Heisenberg kagome AFM remains 
unknown, all of the guidance obtained from VBC or RVB 
constructions suggests a total singlet. Indeed, the nearest- 
neighbor RVB basis has been shown to deliver a semi- 
quantitative account of some of the properties of the system, 
both without^'-^ and with" impurities, and for this reason a 
non-singlet ground state presents a considerable challenge. 
However, by the nature of its structure, a triangular geom- 
etry involving active d^2_y2 and d^^z^-r^ orbitals, the Cu- 
Cu bonds in ZnCu3(OH)6Cl2 are not centrosymmetric. The 
resulting Dzyaloshinskii-Moriya (DM) interactions-^ induce 
quite generally a small admixture of triplet states into a sin- 
glet ground state, and thus have been discussed'^ as a likely 
explanation for the finite response observed at zero tempera- 
ture. The DM interactions in herbertsmithite have been de- 
termined recently by electron spin resonance experiments)^ 
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These measurements show that the dominant component of 
the DM interaction is that perpendicular to the kagome plane, 
which is of order ~ 15 K ~ 0.08J, while the in-plane 
component, D' , although not excluded by symmetry, is much 
smaller (D' ~ 2 K - O.OIJ). 

Clearly, a complete theoretical description for the ground- 
state properties of ZnCu3(OH)6Cl2, including the magnetiza- 
tion response around the impurity sites, must take into account 
the combined effect of both the DM interactions and the non- 
magnetic impurities in the kagome planes. Among the many 
approaches adopted to gain further insight into the kagome 
system, it is well known that the clean, classical kagome 
AFM without DM anisotropy has an extensive ground-state 
degeneracy, which is fully lifted by quantum fluctuations 
only beyond harmonic orderi2^i2i2i The inclusion of DM 
anisotropy perpendicular to the kagome plane selects imme- 
diately the uniform (Q = 0), three-sublattice state,^''" Cepas 
et al. have shown recently^' that in the quantum (s = 1/2) 
kagome AFM, this ordered phase can be stabilized only for 
D/ J > 0.1. This implies that herbertsmithite may well be 
very close to the critical point where quantum fluctuation ef- 
fects promote a disordered ground state. The linear suscep- 
tibility of the clean kagome AFM in the presence of DM in- 
teractions has been analyzed^*^ using a perturbative expansion 
about a short-ranged VBC ground-state scenario, while the ef- 
fect of non-magnetic impurities without DM interactions has 
been studied in Refs. 1,1 L 12.3 1.321 . It has been shown that 
one of the central consequences of a non-magnetic impurity in 
the quantum kagome AFM is a characteristic dimer freezing, 
which takes place around the impurity site due to frustration 
relief in the affected triangle s»iiii^ 

Here we present an extensive exact-diagonalization study 
which includes the effects of both DM anisotropy and non- 
magnetic impurities. Our results can be summarized as fol- 
lows. For D/J< 0.06, the effect of frustration relief around 
the impurit y^'''^ causes the formation of strong, local dimers 
whose magnetization response, a staggered moment directed 
along D X B, is qualitatively similar to that of isolated dimers. 
In this regime we find surprisingly large variations in the 
magnitude of the response on the different induced dimers. 
For sites far away from the impurity, we argue that the re- 
sponse must be uniform and much smaller in magnitude. For 
D/J > 0.1, the system enters the Q = 0, 120°, semiclassi- 
cally ordered state,5iill2i and the effect of impurities becomes 
very short-ranged. However, even here we find that the dimers 
next to the impurity remain strong up to values D/J ^ 1, 
meaning that their spins remain correlated with each other, 
rather than participating in the ordered state, for D/J< 1. 

This article is organized as follows. In Sec.|ll]we introduce 
the model Hamiltonian and discuss some of the symmetry as- 
pects of the clusters considered in our exact-diagonalization 
calculations. The results for the local magnetization response 
are presented and analyzed in Sec. |III] In Sec. |IV| we inves- 
tigate the dominant magnetic correlations on the basis of the 
full correlation matrix (the "natural orbital" method). We con- 
clude in Sec.[V]with a discussion of the implications of our re- 
sults for experimental measurements on ZnCu3(OH)6Cl2. We 
also include two appendices which elaborate on the magneti- 



zation response of a minimal (four-site) cluster (App. |A]l and 
on the natural orbital method (App.lBl). 

II. MODEL 

We consider a spin- 1/2 model on the kagome lattice with a 
single, non-magnetic impurity described by the Hamiltonian 

n = JJ2^' ■ ■ ^ - B • S, (1) 

with periodic boundary conditions. The first term is the 
Heisenberg exchange energy between nearest-neighbor spins 
(ij), the second term represents the DM interactions, and the 
last term is the Zeeman energy of the total spin S = J^i 
in a field B. In what follows, we work in the fixed reference 
frame xyz shown in Fig.[T] The field is taken to be in the xz- 
plane at an angle 6 from the z-axis. The DM vectors on every 
bond are taken to be perpendicular to the kagome plane and 
are fixed by the symmetry of the clean system, namely trans- 
lations and Ce rotations around the hexagon centers. We have 
chosen the site-labelling convention, denoted in Fig.[T]by the 
directionality of the arrows from i to j, such that D^j = De^ 
with D > Q for all bonds. In this study we do not consider 
an in-plane DM component, because in ZnCu3(OH)6Cl2 this 
is known^'* to be much weaker than the out-of -plane compo- 
nent. We note, however, that the in-plane DM problem on the 
kagome lattice is quite significantly different from its out-of- 
plane counterpart studied here.^ 

We will focus on the magnetization response of the sys- 
tem in finite fields, up to the values B lOT (B < J/20) 
probed by experiment. For a general field orientation 9, the 
U(l) spin rotation symmetry around the z-axis is broken and 
thus the total magnetization is not a good quantum number 
Because of the impurity site, the same is true also for the mo- 
mentum. Thus we can treat kagome clusters with a maximum 
of iV = 26 magnetic sites, and here we show primarily the 
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FIG. 1 : (Color online) The kagome lattice of Cu^'*' ions (s = 1 /2), 
containing a single, non-magnetic Zn^^ impurity. Bond arrows de- 
fine site labels i and j in each interaction term D ■ s,; x Sj, while 
D = Dbz with D > for all bonds. Curving arrows denote the Ce 
symmetry about the hexagon centers. 
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(a) D/J=0.05 



(b) D/J=0.1 
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FIG. 2: (Color online) In-plane magnetizations and bond spin cor- 
relation functions for the 26-site kagome cluster with B = J/20 
and — 30°. The lengths of the arrows are proportional to 
((sf)^ + {s^}^)^^^ and the thicknesses of the bonds to (si ■ sj). 
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FIG. 3: (Color online) D-dependence of (a) (sf ), (b) (sf ), (c) {si ■ 
Sj),and(d) (si xsj)^ for the 26-site kagome cluster with B = J/20 
and 9 = 30° . i and j in panels (c) and (d) are nearest-neighbor sites 
only. The curves which correspond to the sites or bonds closest to 
the impurity are designated by (red) filled symbols. 



in the thermodynamic limit, and our results indeed conform 
with this picture. 



results for N — 14, 20, and 26. These clusters are symmetric 
with respect to spatial inversion through the impurity site. For 
a non-degenerate ground state, this sets the constraints 

(Si) = (Sg), (Sa) = (st), (Sg) = (S12), ... (2) 

Furthermore, the clusters with = 14 and 26 are also spa- 
tially symmetric under reflection through the xz-plane which 
passes through the impurity site (Fig.[TJ, whence their Hamil- 
tonian ([U remains invariant under the corresponding mirror 
operation (spatial reflection followed by the time reversal). 
For a non-degenerate ground state, this gives in addition 



This mirror symmetry is also present in the clean, infinite 
kagome lattice. Thus if the vacant site were occupied by a spin 
So, Eq. (|3]l would require (sq) — 0. The remaining symme- 
tries of the clean kagome lattice would then enforce {sf) = 
for all sites i, an argument we will employ below to demon- 
strate that in the disordered phase {D/J < 0.1) there can be 
no staggered magnetic response along the y-axis sufficiently 
far from a vacant site. 

In addition to these "symmetric" clusters, we have also in- 
vestigated clusters with odd numbers of magnetic sites, par- 
ticularly with = 17 and 23, whose ground state has a finite 
moment Sz = 1/2. As shown in Ref. ifssll . this moment does 
not form a bound state around the impurity site, as is the case 
in systems such as the checkerboard lattice,^-' but is delocal- 
ized over the entire cluster. Thus we expect no difference be- 
tween clusters with even or odd numbers A^ of magnetic sites 



III. LOCAL RESPONSE 

Figure |2] shows the local in-plane magnetizations {s^/^) 
(red aiTows) and the local spin coiTelation functions (s^ • Sj) 
(described quantitatively by the thickness of each bond ij) for 
the 26-site cluster with B J/20, 9 = 30°, and for four 
representative values of D/J. The full J?/ J-dependence of 
the in-plane magnetizations, the correlation amplitudes, and 
the local twist amplitudes ((s^ x SjY) is shown for each type 
of site or bond in Fig. [3] These results demonstrate that there 
exist two primary regimes with qualitatively different mag- 
netic response, whose properties we discuss in detail in the 
following subsections. In summary, these are (i) the small- 
D/ J regime (D/J < 0.06), where there is a local, dimer- 
like response for the sites around the impurity, and (ii) the 
large-D/J regime {D/J > 0.1), which is characterized by 
the expected*^' '"'^^ in-plane, 120° magnetic ordering pattern, 
albeit with a special response from the sites right next to 
the vacancy. We also find an apparent intermediate regime 
(0.06 < D/J < 0.1) with a peculiar magnetization pat- 
tern around the impurity, but with no evidence of long-ranged 
magnetic correlations, and we will also characterize this to the 
extent that our calculations allow. 



A. D < 0.06J: local dimer-like response 

We begin by considering the behavior of the nearest- 
neighbor spin correlations, (s^ • s^), which are shown in 
Fig- He) and are also denoted by the bond thicknesses in 
Fig. 121 A low-D regime, where the correlations are controlled 



by the exchange interaction (J) and by the presence of the im- 
purity, can be seen to exist for D/ J < 0.06, although in fact 
the correlation amplitudes begin to deviate from their zero- 
D limit at D/J ~ 0.03. This limit was studied in detail 
m Ref. liltl and is characterized by a large and oscillatory 
variation of the correlation amplitudes as a function of dis- 
tance from the impurity. As shown in Fig. |2a), the strongest 
dimerization is found on the bonds nearest the vacancy, where 
{si ■ Sj) ~ —0.69 is very close to the value of —0.75 ex- 
pected for an isolated dimer. This local "dimer freezing" is a 
key consequence of the frustration relief in the two triangles 
containing the vacancy" '^ and, as we will discuss below, has 
direct implications for the magnetization response of the cor- 
responding sites in the small- _D / J regime. 

The local twist amplitudes ((s^ x Sj)^) shown in Fig.jSjd) 
represent, from Eq. the DM contribution to the total mag- 
netic energy. As expected for the small D/J regime, these 
scale linearly with D/J for all bonds, although with a quite 
inhomogeneous distribution of slopes depending on the dis- 
tance of the bond from the impurity. In particular, the bonds 
closest to the impurity have one of the smallest twist ampli- 
tudes, as might be anticipated from their strong dimerization. 

Turning next to the magnetization response shown in 
Figs. Eta) and [3la)-(b), this develops predominantly in the 
xy-plane, as expected given the easy-plane character of the 
DM anisotropy, and the induced spin components are of order 
i^i'^) ~ 10~^-10~^. By contrast, the out-of-plane response 
is very small ((sf ) < 10^^), and thus will not be consid- 
ered in the following. With the exception of sites strongly 
influenced by the periodic boundary conditions, the magne- 
tization response shows a number of features which are also 
present in the 20-site cluster and can be understood qualita- 
tively on the basis of the strong local dimer formation around 
the impurity. Indeed, the direction of the induced magnetiza- 
tions at the strongly dimerized bonds is qualitatively similar 
to the response expected^"* if they were isolated: a staggered 
polarization in the direction D x B (here the y-axis) and a 
uniform polarization in the direction (D x B) x D (a;-axis). 
At the quantitative level, the actual magnitude of the mag- 
netization response goes beyond the isolated-dimer picture, 
showing the striking feature that the moments induced on the 
four "next-most-strongly" dimerized bonds [bond (3,4) and its 
equivalents in Fig. [T| are larger by a factor of 3-4 than on the 
bonds closest to the impurity. The origin of this result, which 
we probe in detail below, is rather deeper than a simple an- 
ticorrelation between spin correlation amplitude and induced 
moment. It lies in the fact that, while the ground-state wave 
function close to the impurity can be described as an approx- 
imate product state of local dimers, the actual nature of the 
excitations may nevertheless differ significantly from such a 
picture. 

To clarify this point, we isolate the sites 1-6 of Fig. [T] and 
consider their response in the small-!?/ J regime. Figure |4] 
shows the exact magnetization response of this cluster with 
parameters B = J/20, 9 = 30°, and for the same field orien- 
tation and four representative values of D/J shown in Fig. |2] 
The response for D — 0.05 J is very similar to that of the sites 
1-6 in Fig. 12] and in fact the same physics can also be shown 




FIG. 4: (Color online) In-plane magnetization pattern (red arrows) 
and bond strengths (bond thickness in blue) on the isolated six-site 
cluster discussed in Sec. [nil for B = J/20 and 9 = 30°. 



to dominate the response of a four-site cluster (sites 1-4 of 
Fig.[T]i, as detailed in App.lAl In all cases, the magnetizations 
induced on bond (3,4) are larger than those on bond (1,2) by a 
fixed factor of up to 5. The qualitative interpretation (App.lAl) 
of this result lies in the different character of the local exci- 
tations on bonds (1,2) and (3,4): a triplet excitation on bond 
(1,2) is a true eigenstate of the isotropic Hamiltonian, whereas 
a similar triplet excited on the bond (3,4) is not, and may prop- 
agate away. As a result, the magnetization response of bond 
(1,2) is to first order in D/J the same as that of an isolated 
dimer, while that of bond (3,4) will be almost 5 times larger in 
the extreme case of the four-site cluster The same arguments 
apply not only on the six-site cluster, but also to the impure 
kagome lattice, although in this latter case a triplet excitation 
on the bonds next to the impurity does not remain fully local- 
ized because the bonds are not fully "frozen." More generally, 
the qualitative agreement between the response of the small- 
est clusters and the response of the corresponding sites in the 
20 and 26-site calculations shows that, in the regime of small 
D/J, the magnetization profile around the impurity is gov- 
erned almost entirely by local correlations. 

The magnetization response far from the impurity is ex- 
pected to be that of the clean kagome lattice with DM interac- 
tions. As explained above, the staggered magnetization must 
vanish by symmetry, while a uniform component (along the 
X-axis) can survive. This latter scales as D x (D x B) (see 
also Ref. |30]) and thus is very much weaker than the stag- 
gered response in close proximity to the impurity. This re- 
sult is corroborated by exact-diagonalization calculations on a 
clean, 12-site kagome cluster (not shown here), which reveal 
a very weak uniform response ((sf) ~ 10^^). We comment 
that results for larger clean clusters {N =18 and 24) suffer 
from severe finite-size effects due to the presence of inequiv- 
alent loops in these geometries, and thus cannot be used for 
a quantitative determination of the magnetization response far 
from the impurity in the thermodynamic limit. 
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B. D> O.IJ: semiclassical, Q = 0, 120° ordered state 

At large values of D/ J [Fig.|2d)], the magnetization pro- 
files correspond closely to the expecte d^'^"'^^ (5 = 0, semi- 
classical, three-sublattice (120°) state. We note here that in 
this case the magnetic field and the impurity combine to pin 
one of the (infinitely many) degenerate semiclassical states. 
This can be understood physically from the fact that the se- 
lected configuration is the one in which the spin missing at the 
impurity site would be antiparallel to the field, and thus would 
contribute a positive Zeeman energy. However, Fig.|2td) also 
shows that the spins closest to the impurity are slightly tilted 
away from a perfect 120° orientation. In fact these spins do 
not participate in the ordered state in the same way as all the 
other ("bulk") spins do, as can be inferred from the spin corre- 
lation amplitudes in Fig.fSjc). Above D/ J ^ 0.06, and espe- 
cially beyond D/ J ^ 0.6, all the bonds except the two closest 
to the impurity attain very similar strengths. By contrast, the 
dimers next to the impurity (denoted by filled symbols) resist 
this tendency, their spins remaining strongly correlated up to 
surprisingly high values in excess of 13 ~ J. There is no 
such feature on the neighboring bonds, i.e. there is no indica- 
tion of a length scale in this behavior. This is another primary 
consequence of frustration relief on the triangles hosting the 
impurity, and will be further corroborated by our analysis of 
the in-plane correlation data in Sec. lIVI 

One further comment regarding the behavior of the twist 
amplitude at large D/ J is in order here. In the clean kagome 
system, the twist amplitude is the same on every bond, but ac- 
cording to Fig.[3jd), the sites next to the impurity have a twist 
amplitude, ((s x SjY) ~ —0.3 for 13/ J oo, larger than on 
any of the other bonds. This reflects the fact that the relative 
orientation of these two spin pairs lies between the 120° orien- 
tation of the spins on triangles and the 90° orientation favored 
in an isolated dimer, and as such shows the consequence of 
frustration relief for the DM interaction. 



C. Intermediate regime: 0.06 < D/J <0.l 

Our results for the local magnetic properties (Figs.|2]and[3]l 
imply the presence of an intermediate regime, 0.06 ^ D/J < 
0.1, which cannot be classed as being in either of the above 
limits. Here the spins develop a peculiar magnetization pat- 
tern around the impurity site, seemingly indicative of spins 
with no effective interaction which only follow the applied 
field. This is shown for the case D/J — 0.1 in Fig.|2c), but 
the analogous picture for D/J = 0.07 is practically indistin- 
guishable. We have investigated this regime in considerable 
detail as a function of D/J, of the applied field B/J, and of 
the field angle 9, finding essentially perfect linear response to 
the field down to the lowest energy scales. As we will show 
in Sec. |IV] the in-plane magnetic correlation data show no ev- 
idence for long-ranged magnetic order in this regime. 

While we cannot exclude the possibility of some type of 
exotic order in this regime, such as a chiral or a spin-nematic 
phase, we are unable to find any evidence of a local order pa- 
rameter Thus we feel that, pending further studies by other 



approaches, the most likely explanation for the peculiar mag- 
netization patterns we observe lies in finite-size effects. One 
of the reasons that these effects are strong here is the fact that 
a uniform magnetic field is not the appropriate conjugate field 
for the three-sublattice ordered state, and thus competes with 
this type of order close to the critical point where the size of 
the ordered moment is small. 



IV. MAGNETIC CORRELATIONS 

Although the analysis of the local magnetization data in 
Sec. Ill reveals several major, qualitative features of the 
ground state of Eq. ([TJ, a number of issues remain unresolved. 
One is whether there exists any type of magnetic order in the 
window 0.06 < D/J < 0.1, and what is the actual criti- 
cal value of D/J where the 120° phase is established. An- 
other concerns the fact that the magnetization patterns shown 
in Fig. |2] give few indications regarding the consequences of 
finite-size effects in our calculations, which may be quite dif- 
ferent in different parameter regimes, and a method is required 
by which to make this more systematic. Further, the actual 
magnitudes of the moments in the ordered phase cannot be in- 
ferred from the local magnetization data, because the latter do 
not represent the thermodynamic limit. While we will obtain 
these moments, here we will not pursue the observation, from 
our numerical data for the field-dependence of the magnetiza- 
tions at large D (D/J > 1), that a very large field (B/J ~ 1) 
is required to establish their full lengths. 

It is thus necessary to analyze the magnetic correlations in 
the ground state \'^) of Eq. ([T]l as a function of D/J. The 
easy-plane character of the DM anisotropy allows us to focus 
on in-plane correlations, but the breaking of translational in- 
variance by the impurity makes it necessary to consider the 
full set of these, which is contained in an x matrix with 
entries 

a,^{^\s+sj\^). (4) 

We remark here that an analysis based on the connected (or 
cumulant) correlation matrix instead of on C (|4| gives the 
same results because, for the fields considered in this study, 
the local magnetic moments (sf) are negligibly small com- 
pared to the connected portion of the correlations, a result 
which is true for both small and large values of 13/ J. A fur- 
ther consequence of this is that our low-field correlation data 
(B = J/20 in the calculations to follow) are essentially iden- 
tical to the B = Q case. 

For translationally invariant systems, the correlation ma- 
trix can be partially or fully diagonalized by a simple Fourier 
transformation. In the present case, where translational in- 
variance is absent due to the non-magnetic impurity, C must 
be diagonalized numerically. The information contained in 
C is very useful, especially if the ground state contains long- 
ranged magnetic correlations in the a;?;-plane. In this case, it 
is known (Refs. |36,37] and App.lBli that the maximum eigen- 
value Am of C is macroscopically large, i.e. Am oc N. The 
corresponding eigenvector, v,n, then describes the magnetic 
mode of the system with the dominant fluctuations, and thus 
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FIG. 5: (Color online) Eigenvalues Ai of the correlation matrix (di- 
vided by A'^) as a function of D j J for A'^ =14 (circles), 17 (down 
triangles), 20 (squares), 23 (up triangles), and 26 (diamonds). Here 
B = J/20 and 6^ = 30°, but the results for B = are identical. 
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FIG. 6: (Color online) Scaling with system size of the largest eigen- 
value Am (divided by A'^) of the correlation matrix. The solid lines 
correspond to the expected l/\/iV scaling of the leading corrections 
to the thermodynamic limit in the 120° ordered phase. 



is related directly to the spatial dependence of the relevant or- 
der parameter. In the case at hand, the dominant mode is 
not simply a periodic modulation of the spins, but will contain 
valuable information for the non-trivial magnetization profile 
around the impurity. In what follows we exploit this informa- 
tion. 

Figure |5] shows the full eigenvalue spectrum of C for clus- 
ters with N=I4, 17, 20, 23, and 26. The results for B = 
and B = J/20 are indistinguishable. An essential feature of 
Fig.|5]is that, as D j J is increased in the region beyond 0.06, 
the maximum eigenvalue on each of the clusters becomes pro- 
portional to N , and thus very much larger than the remaining 
eigenvalues; we caution that the actual value of this crossover 
cannot be inferred from the data of Fig.|5] where it is evident 
that the curves are still some way from the thermodynamic 



limit, and address this point below. As described in App. |B] 
this means that the system has developed long-ranged in-plane 
magnetic correlations in the regime of large D/J. For small 
D I J , all the eigenvalues depend only weakly onD/ J and are 
closely spaced in magnitude, which is a sign of short-range 
correlations dictated not by D but by J. 

The finite-size scaling of the dominant eigenvalue A„j, nor- 
malized by N, is given in Fig. |6] The extrapolated values 
of this quantity represent the square of the average magnetic 
moment in the thermodynamic limit, where the spins reach 
approximately 80% of their full moment as D/ J oo?'^ For 
all values D/J > 0.1, it is clear that the finite-size correc- 
tions scale as 1/\/N, as expected for a state of broken U(l) 
symmetry.^'' This scaling procedure represents the appropri- 
ate means of deducing the existence of long-ranged magnetic 
order, by continuing the curves of Fig.|5]to the infinite-system 
limit. However, this powerful method shows no indication 
of such order in the regime 0.06 < D/J < 0.1, specifying 
that the transition to the semiclassical state should be taken as 
D/J- 0.1. 

The magnetization profile corresponding to the dominant 
eigenmode Vm also contains important information, which is 
shown in Fig. |2]for the four representative values of D/J and 
represented by two-dimensional arrows whose components 
are the real and imaginary parts of v^- At D/J < 0.06, 
there is no dominant mode as is the case at large D, but the 
strongest mode shown in Fig. |2la) corresponds nevertheless 
to the pattern of strong spin correlations around the impurity 
[Fig.|2fa)]: the correlations in this mode are confined to the 
strong bonds next to the impurity, where the spins are almost 
antiparallel. The strength of these local correlations is gov- 
erned by J, which is the reason why X„i remains essentially 
D-independent for D/J< 0.06 in Fig. |5] We emphasize 
again that the profile shown in Fig.|2la) does not represent the 
actual magnetization response - this is shown in Fig. - 
but rather the dominant fluctuation mode. 

The situation changes dramatically at D/J > 0.1 
[Figs.|3c)-(d)], where the system develops long-ranged order 
with the majority of spins participating in the Q = 0, semi- 
classical 120° state. The data show clearly that the crossover 
from the dimer-Iike regime at small D/J to the ordered phase 
at high D/J is smooth, and there are no indications for 
any other type of magnetic order in the intermediate regime 
0.06 < D/J < 0.1. This suggests once again that the pecu- 
liar magnetization pattern found in Fig.|2|b) is most likely to 
be a fluctuation effect which is emphasized close to the critical 
point. 

Indeed, our analysis of the magnetic correlations allows us 
to develop a qualitative understanding of the situation in this 
regime. The kagome AFM is characterized by a very high 
density of nearly-degenerate states in the ground manifold. 
The effect of a DM interaction term D is to favor certain 
states and penalize others, but this process is extremely sen- 
sitive to the shape and size of the cluster used in the calcula- 
tions. Some of our results, such as the discrepancy between 
Figs.|2c) and|7Jc), are evidence of significant finite-size ef- 
fects. These effects are expected to, and indeed in Figs. |5] 
and|6]observed to, increase on approaching the critical regime 
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around _D/ J ~ 0.1. With such small separations between sets 
of energy levels, the effect of a magnetic field applied in this 
regime can be dramatic: certainly some of the differences be- 
tween Figs. |2] and |2] may be ascribed to the finite field in the 
former While it remains impossible to exclude more exotic 
types of order for 0.06 < D/J < 0.1, whose fingerprints 
might be difficult or even impossible to discern in data ob- 
tained from the cluster sizes accessible by exact diagonaliza- 
tion, these considerations make systematic our statement that 
the most probable explanation for the tendency of the small- 
D/ J and large-Z)/ J data to suggest the presence of an inter- 
mediate regime remains in the realm of finite-size effects. 

The magnetic correlation data also quantify the extent to 
which the four sites closest to the impurity continue to resist 
the effect of D (Fig.|7]). As already inferred from the behavior 
of the correlation amplitudes (Fig.|2]i, these spins do not adopt 
a relative 120° orientation, but remain almost antiparallel to 
each other up to DM interactions D ^ J. This direct con- 
sequence of frustration relief in the two doped triangles also 
suggests that the actual orientation of the ordered component 
of these spins will continue to reflect the behavior of a dimer, 
and thus will be dictated by D and B as well as by the cou- 
pling to the other spins. Only at surprisingly high values of 
D/J does a crossover occur to predominant correlation with 
the other spins of the system, and hence to participation in the 
120° ordered state. 

This special response of the four sites closest to the impu- 
rity is demonstrated once again in Fig. |8] which shows the 
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FIG. 7: (Color online) Magnetization profile for the 26-site kagome 
cluster in zero field, corresponding to the eigenvector Vm of C with 
the largest eigenvalue Am . The in-plane moments are given by the 
real and imaginary parts of Vm- Note that this mode is unique up to 
a global U(l) rotation which is related to the arbitrary phase of Dm. 



FIG. 8: (Color online) In-plane spin magnitudes and orientations 
(represented respectively by \/Xm\v,n{i)\ and by the phase 8i of 
Vm{i)) as a function of D/J, for clusters with A'^ = 14 (circles), 
20 (squares), and 26 (triangles). The data corresponding to the four 
sites closest to the impurity are highlighted by filled symbols. 



at site i is represented by ^/)wn\v„i{i)\, while the correspond- 
ing angle is determined from the phase representation of the 
real and imaginary parts of The primary feature of 

Fig.[8]is that, while both the magnitude and the orientation for 
the majority of sites (open symbols) largely merge together at 
D/J> 0.06, the same quantities for the four sites closest to 
the impurity (filled symbols) remain different up to Z? ^ J. 
This feature, present in all three clusters, reinforces the picture 
of the special behavior of the frustration-reheved spins. 



V. SUMMARY AND DISCUSSION 

We have presented an extensive exact-diagonalization study 
of the combined effects of non-magnetic impurities and 
Dzyaloshinskii-Moriya (DM) interactions in the s = 1/2 
kagome antiferromagnet (AFM). Although the clean and 
purely Heisenberg kagome AFM remains an enigma, being 
apparently non-magnetic and possibly dimer-based, we have 
shown how DM interactions immediately induce a magnetic 
response, and impurities immediately nucleate very strong 
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and somewhat extended patterns of dimerization. We have 
found that the system undergoes a phase transition from this 
magnetically disordered state to the Q = 0, three-sublattice 
(120°) ordered state when the interaction ratio between the 
DM and APM terms is D/J ^ 0.1. This conclusion is in 
agreement with the results obtained by Cepas et al. for the 
clean kagome lattice. 

For D/ J < 0.06, there is strong dimerization of the spins 
next to each impurity site as a consequence of frustration re- 
lief in these two triangles, as shown by Dommange et alM for 
the D = case. An oscillating pattern of weakly and strongly 
dimerized bonds then extends some distance from the impu- 
rity site. In the presence of DM interactions, the applied field 
induces a magnetic response on top of this dimerization pat- 
tern, which on the strongly dimerized bonds shows the same 
qualitative behavior as an isolated dimer, namely a dominant 
staggered component directed along D x B. The induced 
magnetic moments on the sites closest to the vacancy are small 
as a result of their strong dimerization, and we have shown 
why the magnetization on the next-nearest strong bonds is ap- 
proximately four times larger than this. The maximum in- 
duced magnetization is thought to occur on these sites or on 
the next (i.e. next-next-nearest) "ring" of more strongly dimer- 
ized bonds beyond them. ' ' Still further from the impurity, the 
local moment decreases again, towards the weak and uniform 
response (by which is meant no staggered component) of the 
clean kagome system with DM interactions. Because this is 
proportional to D x (D x B), it is much smaller than the 
staggered response around the impurity sites. 

The ordered phase established at D/J > 0.1 shows all the 
hallmarks expected for semiclassical magnetic order How- 
ever, we found that the spins closest to the impurity remain 
strongly correlated with each other up to surprisingly large 
values, D ^ J, before they begin to participate fully in the 
ordered state. Thus a substantial portion of the actual mag- 
netization of these spins remains similar to that of a dimer, 
and is set by D and B. Despite the fact that the obvious 
features of these two phases end at separate values of D/J 
which leave open the possibility of an intermediate phase, be- 
tween these limits we have been unable to find any definite 
evidence for another type of physical behavior Thus we are 
forced to the tentative conclusion that there is only one true 
phase transition in the system, and the features we have ob- 
served for 0.06 < D/J < 0.1 are most likely to be artifacts 
resulting from the low effective energy scales in this regime 
and from the small system sizes to which our calculations are 
constrained. 

We conclude by discussing the implications of our re- 
sults for the interpretation of experiments performed on 
ZnCu3(OH)6Cl2. The most accurate probe of the local mag- 
netic response around doped impurities is NMR, and so we 
focus on the ^^O NMR measurements of Olariu et al. Be- 
fore comparing the experimental results to our calculations, 
which were performed for a single vacancy in small kagome 
clusters, it is necessary to consider the impurity concentra- 
tions in the experimental samples. Because these are on the 
order of 5% for the cleaner samples used in the measurements 
made to date,'^ the average linear spacing between impurities 



is then 4-5 Cu-Cu bond spacings: in fact this is rather simi- 
lar to the lengthscales probed in our N = 20 and 26 clusters. 
From our results, it is safe to say that there are no Cu spins 
in the real material which can be considered as being far from 
any impurity sites. 

Our conclusions show that the situation D/J > O.lmaybe 
safely excluded for herbertsmithite, because a magnetically 
ordered phase of any type is inconsistent not only with the 
NMR experiments but also with muon spectroscopy and in- 
elastic neutron scattering experiments,'^ '^ where no such or- 
der has been found for any temperatures down to 50 mK. By 
contrast, for DM interactions D/J < 0.1, one would expect to 
observe a broad NMR signal as a consequence of the consider- 
able number of inequivalent Cu sites produced by the partially 
frozen patterns of dimerization and induced local moments 
around the randomly distributed impurities. This is compati- 
ble with the shape of the (M) line measured in Ref. [18], while 
the (D) line, which probes only one type of Cu site, would in- 
deed be expected to be sharp. Further, our calculations quan- 
tify the way in which the DM interaction removes the apparent 
singlet nature of the D — ground state, leading to the finite 
lineshift observed in experiment in the limit T ^ 0. 

To be more specific about the actual strength of D/J in her- 
bertsmithite within the disordered regime, we appeal to our 
results (Sec. Ill) for the relative local moments on the sites 
around an impurity. On the qualitative assumption that the 
^^O line shift should scale with the sum of the moments on 
the two neighboring Cu sites, the 1:4 moment ratio found for 
the nearest two Cu sites in the low-_D regime (D/ J < 0.06) 
would be expected to result in a 1:5 ratio of line shifts. Thus 
the observed 1:2 ratic''^ suggests quite strongly that herbert- 
smithite falls in the regime 0.06 D / J < 0.1 while the exact 
nature of the ground state in this region remains to be con- 
firmed, it is clear that it is non-magnetic, and that the moment 
ratio between the two types of Cu sites closest to an impurity 
is significantly smaller than at low D/J [Fig. |8jb)]. Such a 
value for D/J is fully consistent with the result D ~ 0.08 J 
obtained by electron spin resonance in Ref. i24ll . 

Finally, we have shown in the disordered regime {D/J< 
0.1) that sites close to the impurities exhibit strong features 
characteristic of the response of isolated dimers, with those 
right next to an impurity (the (D) line) reflecting a very weak 
induced magnetization. This behavior could explain the sup- 
pressed local susceptibility at T ^ and the enhanced spin- 
relaxation times Ti found'^ for these sites. However, to go 
beyond this qualitative level of agreement, it seems neces- 
sary to have access to single-crystal data, which would pro- 
vide a detailed understanding of the hyperfine interactions in 
ZnCu3(OH)6Cl2. These are required to investigate the depen- 
dence of the magnetic response on the field orientation and 
the line-shift contributions from each Cu site as a function of 
its induced magnetization, information which could be inter- 
preted directly within our calculational approach. 
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FIG. 9: (color online) Typical in-plane magnetization pattern (red ar- 
rows) and spin correlation functions (represented in blue by the bond 
thickness) for the minimal four-site cluster in the small-D regime. 
Here D = B ^ J/20, 9 = 30°, and the orientations of the DM 
interactions are the same as for sites 1-4 in Fig.[T] 
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APPENDIX A: LINEAR RESPONSE OF THE FOUR-SITE 
CLUSTER 

Here we present the details of the magnetic response of the 
four-site cluster shown in Fig.|9] and demonstrate why this can 
be considered as the minimal cluster illustrating the physics of 
the low-D regime. The Hamiltonian is 

n^ni + H2 + v = Ho + v, (Ai) 

in which 

Til = J(si • S2 + S3 • S4) - B • S, (A2) 

H2 = J(S2 • S3 + S2 • S4), (A3) 

V = D • (s2 xsi + S2 XS3 + S3XS4 + S4XS2), (A4) 

S = J2i is the total spin, D = De^, and B lies in the 
x2-plane, as shown in Fig. [T] Our purpose is to examine the 
magnetization response of this cluster for small fields and to 
explain why this captures the generic response of the corre- 
sponding sites around the impurity in Fig.|2] We find numeri- 
cally that to leading order inD/J there is a staggered response 
along the y-axis given by 

(si) = -(s2) = iDxB (A5) 

(S3) = -(s4) = yDxB, (A6) 

i.e. the magnetization response at bond (1,2) is precisely that 
expected for an isolated dimer while, surprisingly, the re- 
sponse at the bond (3,4) is approximately five times greater. 
We provide a qualitative interpretation of this result based on 
the character of the excitations of Ho . 

Using V as a perturbation, the ground state is given to lead- 
ing order by j?/') = (1 + 7^V)|V'o). where |i/'o) is the ground 
state of TYq, TZ = is the corresponding resolvent oper- 

ator with Q — 1~ |V'o)('0o|i and — —3 J/2. The mean 



value of any operator A is then given by 

(A) - (Vo|A7^V|Vo> +H.C., (A7) 

where we have assumed that {■4>q\A\iI)q) — 0. Because this 
relation gives the linear response, it may be used to obtain the 
individual contribution of any of the DM terms contained in 
V to the magnetization of any spin site. 

For each of the bonds (1,2) and (3,4), it is convenient to use 
the singlet and triplet basis states \s) (| ti) - | iT))/\/2, 

\ti) = I TT), |io) = (I Ti) + I iT))/V2, and = | [[). 
For B < J/2, the ground state \iPq) of Hq is the product state 
of the two singlets, 

Ho) = \s) 12(E) \s)u- (A8) 
Also applicable are the relations^ 

H2|s)l2® |s)34-0 (A9) 

H2|i™)i2<8 |s)34 = 0, Vto = 0,±1, (AlO) 

which specify that the state |t„i)i2 'S) {3)34 is also an eigenstate 
of Ho- Thus the local excitation of a triplet on bond (1,2), for 
example by acting with the term D • S2 x si (or the operators 
A = 5^ 2) on IV'o), gives 

7^|tm)l2 «) |s)34 = —^—^\tm}l2 » |s)34, (All) 

1 + niB 

and the triplet remains localized. 

By constrast, a triplet on bond (3,4) will not remain local- 
ized, because |s)i2 (8) |tm)34 is not an eigenstate ofHo- With 
these considerations and by using Eq. (|A7| i, it is easy to show 
that of all the DM terms contained in V, only D • S2 x si con- 
tributes to (si_2), and this explains in turn why the response 
at bond (1,2) is that of an isolated dimer. The staggered mag- 
netization at bond (3,4), on the other hand, is driven by all 
the terms of V except D • S2 x Si, which explains why this is 
different from the response of bond (1,2). More generally, we 
see that although the ground state of Ho is a product of two 
dimers, the different character of the local excitations leads 
to quantitatively very different responses. This completes the 
qualitative interpretation, summarized in Sec. [Till for the dif- 
ferent magnetization responses of the sites around the impu- 
rity in Fig.|2] 



APPENDIX B: INTERPRETATION OF MAGNETIC 
CORRELATIONS 

Here we provide a more detailed presentation of the natural 
orbital method used in Sec.|IV]for the study of in-plane mag- 
netic coiTelations. We consider the ground state \'^) of the 
Hamiltonian given in Eq. ([T]) and express the in-plane mag- 
netic coiTelations by the matrix C of Eq. (01). Let us denote by 
{^a;Va} the set of eigenvalues and normalized eigenvectors 
of C, i.e. C ■ Va = XaVa- In addition to being Hermitian, the 
correlation matrix C is also positive semi-definite, because for 
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any normalized vector one has 1)1^ • C • = |lA^i,|^')|p > 0, 
where 

M„ = ^t;(z)s- (Bl) 

i 

defines a macroscopic magnetic mode. Thus > for all 
a. Further, because s+s^ = 1/2 + for s — 1/2, TrC — 
Aq = N/2 + (5^). We note also that the eigenvalues Aq 
give the fluctuations of the matrix M.a (specified by Eq. (IB lb 
with V = Va), because 

{^MlMl>\^) = XaSab- (B2) 

In order to discuss magnetic order, we consider the thermo- 
dynamic limit and address the question of what is required to 
ensure long-ranged in-plane magnetic correlations, i.e. 

lim C,j ^ 0. (B3) 

To this end, it is convenient to express C in terms of its spec- 
tral decomposition, 

C,, =^Aa«a(zK(j). (B4) 

a 

Equations ( IB3l l and (IB4l l mean that there should exist at least 
one eigenstate v„i with a nonzero amplitude v,n{i) for all i. 
However, v is a quantity normalized to one by an overall fac- 
tor of 1/ ^/N, which from Eq. ( IB3l l means that the correspond- 
ing eigenvalue is macroscopically large, A^ cx N. If only one 
(the maximum) eigenvalue has this property, it is safe to re- 
place Eq. ( lB4b with 

Cy ~ XrnVmii)v^{j), (B5) 



because the neglected terms are negligibly small for iV — > oo. 
This special, separable form of the correlation matrix in an or- 
dered state allows one to identify the relevant local order pa- 
rameter. Indeed, in an explicitly symmetry-broken (coherent) 
state we may replace (s^sj) by {s^){sj), and hence make 
the identification 

(4) = (s?> + «(s.r> = Vx^v^it). (B6) 

We remark here that these in-plane spin components are fixed 
only up to a global U(l) rotation, due to the fact that the eigen- 
state Vjn is only specified up to a global phase. The corre- 
sponding macroscopic order parameter is given by Aim, be- 
cause [using Eq. ( IB2l l1 

mMlMml'i') = A™ (X TV, (B7) 

and thus the response to a conjugate field that couples directly 
to Aim diverges in the thermodynamic limit. 

Finally we should mention the analogy of the discussion 
presented in this appendix to Bose-Einstein condensation, and 
in particular to the case of cold-atom systems confined in 
a harmonic trap (see for example Ref. |41]). The correla- 
tion matrix in such a system is the one-body density matrix, 
/9i(x, x'), which measures the coherence between different 
parts of the system. The eigenstates of pi are termed "the nat- 
ural orbitals", and the eigenvalues give the relative occupation 
probability of these orbitals. As in our situation, off-diagonal, 
long-range order is signaled by the fact that one of the eigen- 
values of pi becomes macroscopically large. ^^'^^ The domi- 
nant eigenstate Vm of a magnetic system is thus analogous to 
the condensate wavefunction in a superfluid. 



' G. Misguich and C. Lhuillier, in Frustrated Spin Systems, edited 
by H. T. Diep (World Scientific, Singapore, 2005). 

^ P. W. Anderson, Mater. Res. Bull. 8, 153 (1973). 

^ J. B. Marston and C. Zeng, J. Appl. Phys. 69, 5962 (1992); A. 
V. Syromyatnikov and S. V. Maleyev, Phys. Rev. B 66, 132408 
(2002); P Nikolic and T. Senthil, Phys. Rev. B 68, 214415 (2003); 
R. Budnik and A. Auerbach, Phys. Rev. Lett. 93,187205 (2004); 
R. R. P Singh and D. Huse, Phys. Rev. B 76, 180407(R) (2007). 

" S. Sachdev, Phys. Rev. B 45, 12377 (1992); F. Wang and A. Vish- 
wanath, Phys. Rev. B 74, 174423 (2006). 

^ R Mila, Phys. Rev. Lett. 81, 2356 (1998). 

^ H. C. Jiang, Z. Y. Weng, and D. N. Sheng, Phys. Rev. Lett. 101, 
117203 (2008). 

M. B. Hastings, Phys. Rev. B 63, 014413 (2000). 
^ Y. Ran, M. Hermele, P A. Lee, and X.-G. Wen, Phys. Rev. Lett. 

98, 117205 (2007); M. Hermele, Y. Ran, P A. Lee, and X.-G. 

Wen, Phys. Rev. B 77, 224413 (2008). 
' M. Elhajal, B. Canals, and C. Lacroix, Phys. Rev. B 66, 014422 

(2002). 

R. Ballou, B. Canals, M. Elhajal, C. Lacroix, and A. S. Wills, J. 
Magn. Magn. Mater. 262, 465 (2003). 
" S. Dommange, M. Mambrini, B. Normand, and F. Mila, Phys. 



Rev. B 68, 224416 (2003). 

A. Lauchli, S. Dommange, B. Normand, and F. Mila, Phys. Rev. 
B 76, 144113 (2007). 

M. P Shores, E. A. Nytko, B. M. Barlett, and D. G. Nocera, J. 
Am. Chem. Soc. 127, 13462 (2005). 

M. Rigol and R. R. P Singh, Phys. Rev Lett. 98, 207204 (2007). 

G. Misguich and P Sindzingre, Eur. Phys. J. B 59, 305 (2007). 

P. Mendels, F. Bert, M. A. de Vries, A. Olariu, A. Harrison, F. 

Due, J.-C. Trombe, J. S. Lord, A. Amato, and C. Baines, Phys. 

Rev. Lett. 98, 077204 (2007). 
" J. S. Helton, K. Matan, M. P Shores, E. A. Nytko, B. M. Bartlett, 

Y. Yoshida, Y. Takano, A. Suslov, Y. Qiu, J.-H. Chung, D. G. 

Nocera, and Y. S. Lee, Phys. Rev. Lett. 98, 107204 (2007). 

A. Olariu, P Mendels, F. Bert, F. Due, J. C. Trombe, M. A. de 

Vries, and A. Harrison, Phys. Rev Lett. 100, 087202 (2008). 
" S.-H. Lee, H. Kikuchi, Y. Qiu, B. Lake, Q. Huang, K. Habicht, 

and K. Kiefer, Nat. Mater. 6, 853 (2007). 
^° M. A. de Vries, K. V. Kamenev, W. A. Kockelmann, J. Sanchez- 

Benitez, and A. Harrison, Phys. Rev Lett. 100, 157205 (2008). 
^' F. Bert, S. Nakamae, F. Ladieu, D. L'Hote, P Bonville, F. Due, 

J.-C. Trombe, and P Mendels, Phys. Rev B 76, 132411 (2007). 

M. Mambrini and F Mila, Eur. Phys. J. B 17, 651 (2000). 



11 



I. E. Dzyaloshinskii, J. Phys. Chem. Solids, 4, 241 (1958); T. 
Moriya, Phys. Rev. Lett. 4, 228 (1960). 

A. Zorko, S. Nellutla, J. van Tol, L. C. Brunei, F. Bert, F. Due, 

J.-C. Trombe, M. A. de Vries, A. Harrison, and P. Mendels, Phys. 

Rev. Lett. 101, 026405 (2008). 

R. Moessner, Can. J. Phys. 79, 1283 (2001). 

A. B. Harris, C. Kallin, and A. J. Berlinsky, Phys. Rev. B 45, 2899 

(1992). 

" A. Chubukov, Phys. Rev. Lett. 69, 832 (1992). 

C. L. Henley and E. R Chan, J. Mag. Mag. Mater. 140-144, 1693- 
1694 (1995). 

O. Cepas, C. M. Fong, P W. Leung, and C. Lhuillier, Phys. Rev. 
B 78, 140405 (R) (2008). 

M. Tovar, K. S. Raman, and K. Shtengel, Phys. Rev. B 79, 024405 
(2009). 

K. Gregor and O. I. Motrunich, Phys. Rev. B 77, 184423 (2008). 
M. J. Rosenberg and R. Chitra, Phys. Rev. B 78, 132406 (2008). 

D. Poilblanc, A. Lauchli, M. Mambrini, and F. Mila, Phys. Rev. B 
73, 100403 (R) (2006). 

S. Miyahara, J.-B. Fouet, S. R. Manmana, R. M. Noack, H. 
Mayaffre, I. Sheikin, C. Berthier, and F. Mila, Phys. Rev. B 75, 



184402 (2007). 

S. Capponi, A. Lauchli, and M. Mambrini, Phys. Rev. B 70, 
104424 (2004). 

0. Penrose and L. Onsager, Phys. Rev. 104, 576 (1956). 
" C. N. Yang, Rev. Mod. Phys. 34, 694 (1962). 

We note here that the the extrapolated values of the (unrescaled) 
squared moments reported for the clean kagome system in 
Ref. i29ll are smaller than the true values by a factor of two. The 
origin of this factor lies in the choice of the a;-components of the 
spins for studying the correlations: The dominant eigenvalue of 
these correlations is twofold degenerate and corresponds to two 
modes with the same magnetization pattern in the s-direction but 
with different chirality. By contrast, the correlation matrix C de- 
fined in Eq. (HJl is able to select the state with the right chirality, 

1. e. the one favored by the DM interactions. 

H. Neuberger and T. Ziman, Phys. Rev. B 39, 2608 (1989); D. S. 
Fisher, ibid. 39, 11783 (1989). 

S. Miyahara and K. Ueda, Phys. Rev. Lett. 82, 3701 (1999). 
J. L. DuBois and H. R. Clyde, Phys. Rev. A 63, 023602 (2001). 



